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Abstract 

Filters whose porosity decreases with depth are often more efficient at removing so¬ 
lute from a fluid than filters with a uniform porosity. We investigate this phenomenon 
via an extension of homogenization theory that accounts for a macroscale variation in 
microstructure. In the first stage of the paper, we homogenize the problems of flow 
through a filter with a near-periodic microstructure and of solute transport due to ad- 
vection, diffusion, and filter adsorption. In the second stage, we use the computationally 
efficient homogenized equations to investigate and quantify why porosity gradients can 
improve filter efficiency. We find that a porosity gradient has a much larger effect on 
the uniformity of adsorption than it does on the total adsorption. This allows us to 
understand how a decreasing porosity can lead to a greater filter efficiency, by lowering 
the risk of localized blocking while maintaining the rate of total contaminant removal. 


1 Introduction 


Membrane separation is a vast industry 

treatment [l|,[^, biopharmaceuticals [1,1^, and food processing _ _ _ 

crucial to remove waste and excess water from the blood in kidney dialysis , and yeast and 
bacteria in beer production [^. Despite the diverse industrial applications, an overarching 
goal of membrane design is to maximize the product yield. In applications where the solvent 
is the desired product, such as in water treatment, this is accomplished by maximizing both 
particle removal from the fluid suspension and filter lifespan. Mathematical modelling can 
offer key insight into the filtration process and operating conditions, and thus provide a 
cost-effective way to optimize filter design. 

A common form of membrane separation is depth filtration, in which small particles 
(contaminants) are trapped within, and not just on the surface of, the porous filter material. 
Such filters often capture the majority of the contaminants in the initial portion of the 
filter while leaving the latter portions relatively unused, leading to premature clogging and 
reduced filtration efficiency j^. In such cases, filters whose porosity decreases with depth, or 
porosity-graded filters, which are an example of a functionally graded material, can improve 


with a wide range of applications, including water 

. For example, filters are 


filtration efficiency, and so are often used experimentally lOKlJ]- Their increased efficiency 


can be qualitatively attributed to a decrease in porosity compensating for a reduction of 
contaminant concentration with depth, which occurs due to prior filtering. However, the 
mechanism behind this observation is not fully understood. Experiments are costly and 
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moreover it is very difficult to observe the particle trapping within a filter during the filtration 
process directly. Instead, deductions are made only after dissecting the porous medium once 
filtration has ceased. For these reasons, optimizing a porosity distribution via a systematic 
series of experiments is impractical. 

Mathematical and computational methods to model the filtration process are very useful 
for investigating this issue at a fraction of the full experimental cost. A summary table of 
Computational Fluid Dynamics (CFD) models for filtration is given in 1^. In the same 
paper, the authors use scanning electron microscopes to obtain a full description of a given 
membrane microstructure, and implement a full CFD model of particles moving within the 
membrane. Whilst a full CFD simulation gives excellent insight into how an individual 
particle is trapped, there are large computational costs associated with keeping track of all 
the particles within a complicated pore structure. Additionally, as in the example given 
above, each membrane must be scanned to accurately represent its pore structure. Although 
filtration occurs on the scale of the particle or pore size, it is generally the overall macroscale 
behaviour, such as the total mass of particles removed, which is the main concern in filter 
design. 

An alternative approach to CFD for complex heterogeneous materials, such as filters, is 
to consider upscaling methods. These reduce the complexity of an equation by averaging 
any microscale variation while retaining the important macroscale variation. Mathematical 
homogenization via the method of multiple scales is an asymptotic technique often used for 
this upscaling. Traditionally, for this technique to be appropriate to use, there must be a 
periodic microstructure and the ratio between the length of the periodic cell and the length of 
the important macroscale variation must be small. In fluid mechanics, where the underlying 
flow equations can be difficult to solve, this procedure is used to determine the macroscale 
behaviour of fluids in complicated geometries. For example, in saturated single-phase Stokes 
flow past a periodic array of inert obstacles, this homogenization procedure leads to Darcy’s 
Law 16|, and, in non-saturated media, homogenization techniques can be used to derive 
Richards’ equation |l7l |. 

Whilst the majority of homogenization procedures are carried out in domains whose mi¬ 
crostructure is fully periodic, and thus would only be applicable to filters with a constant 
macroscale porosity, extensions to non-periodic domains have been explored. In 1^, the 


steady problem of nutrient uptake past randomly placed point sinks is considered in one 
spatial dimension. As the governing equations can be solved if the locations of the sinks are 
known, significant analytic progress is made into investigating the macroscale effect of differ¬ 
ent random distributions. A comprehensive extension to near-periodic domains is developed 
in [l^, where the authors use a general curvilinear coordinate transform to homogenize the 
electric potential within a beating heart, mapping the near-periodic microscale to a periodic 
domain. Similar approaches are used to develop homogenized equations in 20|, where the 


authors consider saturated flow in a poroelastic solid with surface growth of the solid phase, 
which is allowed to have a macroscopic variation, and in 21[, where the authors consider 
blood and drug exchange across capillary walls in malignant tumours. 

Although these extensions to non-periodic domains are certainly significant advancements 
to the field, their generality means that, usually, the cell problem must be solved at e very 
point in the macroscale, reducing the computational efficiency. This issue is bypassed in [22l |. 
where the multiple scales homogenization method is extended to quasi-periodic structures 
where the microstructure is allowed to vary slowly. In particular, the authors consider the 
problem of diffusion in a porous medium whose solid phase is modelled as an array of inert 
spherical inclusions (obstacles). The radii of these obstacles are allowed to vary spatially 
over the macroscale length. By imposing a specific one-parameter form on the obstacles, 
an explicit formulation of the macroscale equation is obtained in terms of the cell-averaged 
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porosity, which varies in the macroscale. A notable aspect of this analysis is that a porosity 
variation induces a macroscale particle advection in the direction of decreasing porosity. 

H, 


In this paper we employ the homogenization method developed in 


to investigate 


porosity-graded filters. We model the filter material as a collection of spherical obstructions, 
whose radii vary spatially over a long lengthscale, past which a fluid with suspended contami¬ 
nant particles flows due to an applied transmembrane pressure. The particles are transported 
via advection and diffusion, and can be trapped on the filter microstructure. In general, the 
accumulation of contaminant particles via trapping modifies the filter geometry, and this 
process eventually leads to pore blockage within the filter. Since this effect occurs within 
and not just on the inlet surface of the filter, the blockage is difficult to remove. Dealing 
with a blockage causes a significant reduction in efficiency and an increase in maintenance 
cost, as the filtration procedure must be halted to perform a back-flow procedure to dislodge 
the blockage or, more drastically, a filter replacement. 

The underlying trapping or adsorption mechanism of a particle to the filter structure 
will, in principle, depend on the solute solubility and the filter material adsorption. Many 
existing models for trapping in filters use constitutive macroscale laws to determine how fluid 
volume throughput varies in time, which are fitted with experimental data 2^. In recent 


work by 2^, the authors perform simulations of a microscale model for the blocking of a 


network of pores formulated from a more fundamental mechanistic level, and their results are 
consistent with the set of constitutive laws. However, the model in 2^ does not account for 
the flow profile and is computationally expensive as it keeps track of each pore. Incorporating 
pore blockage and the subsequent filter failure is mathematically challenging to include in 
homogenized models, since pore growth leads to difficult free-boundary problems. 

The main objective of this work is to understand and quantify how porosity-graded filters 
improve filter efficiency. To this end, starting from the flow and particle transport problems in 
the complicated domain described by a given microstructure, we employ the homogenization 
method developed in 19|, [2^ to systematically determine an effective macroscale equation. 
Our model takes into account the effects of the microstructure on the fluid flow, as well as the 
transport by diffusion and advection of the contaminant particles. The trapping of particles 
by the filter microstructure is also included as an adsorption process. As in [2^, we consider 
a solid structure composed of spherical inclusions whose radii vary slowly in the macroscale. 
However, now the inclusions are partially adsorbing sinks for the particle transport and inert 
obstacles to the flow. In order to focus on the effect of variations in porosity within a filter, 
in this work we assume that the contaminant particles are small and in dilute suspension 
within the fluid, and thus their trapping has a negligible effect on the obstacle size. This 
means that our model will not be able to explicitly capture pore blocking. As a result of this 
simplification, the flow problem decouples from the particle transport problem (but not vice 
versa), and we do not have a free-boundary problem. As we consider a dilute suspension, 
we follow the lead of 2^ and impose that particle adsorption is linearly dependent on the 
number of particles. 

This paper is divided into two stages. The first stage consists of the presentation of 
the full problem in ^ and the homogenization of the flow and particle transport problems 
to systematically obtain effective equations on the macroscale in ^ This allows us to 
investigate the effect of a varying porosity on hltration in the second stage of the paper, 
where the steady-state effective equations are applied to a filter whose porosity varies in one 
direction only and the flow is uni-directional, in the same direction as the porosity variation. 
We consider a general operating regime where the filter is placed between two reservoirs. We 
solve the derived system numerically and analyse our results to investigate how a variation in 
hlter porosity affects particle trapping in steady state operation. We introduce two metrics 
to quantify the effectiveness of a filter and, through these, we are able to determine how 
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a porosity gradient can improve filter efficiency. Additionally, we perform an asymptotic 
analysis on the derived system, exploiting the small porosity gradient (a common feature 
in many industrial operating regimes) to determine very accurate analytic expressions at 
a much reduced computational cost. These tasks are all carried out in 0 We conclude 
in ^ with an overview and discussion of the results from this paper, ending with relevant 
extensions based on this work. 

2 Model description 

We consider the transport of particles via advection and diffusion through a porosity-graded 
filter. The filter is modelled as a collection of solid obstacles, to which particles can adsorb. 
We describe the particle locations in terms of the concentration c(®, t), where x is the spatial 
vector coordinate and t is time. 

The concentration field is defined within the fluid phase of the domain, where Qf C is 
the fluid phase and C is the entire domain (which we refer to as the ‘porous medium’), 
with d = 2 or 3 the number of spatial dimensions. We define the solid phase of the domain 
as Qg C noting that Qg = ^\^f- The solid phase is modelled by a collection of fixed 
non-overlapping d-dimensional balls, whose centres are located on a square or cubic lattice at 
a distance dl apart, where d is a (small) dimensionless parameter and I is the characteristic 
filter length. We allow the radii of the balls to vary in space, and a ball with centre at x has 
radius 5lR{x). 

A consequence of the assumption that the contaminants are small particles which are 
in dilute suspension is that particle-particle interactions are negligible, and the dominant 
interaction is the adsorption of particles on the obstacles. Thus, the particle concentra¬ 
tion is governed by the standard advection-diffusion equation with a reactive (or partially 
adsorbing) Robin boundary condition, which represents a linear adsorption rate: 

^ = V • (^DVc — udj , X £ Qf, (2.1a) 

—jc = n‘(^DVc — udj, x£dVlf. (2.1b) 

Here V refers to the nabla operator with respect to H > 0 is the constant diffusion 
coefficient, d^f is the interface between the solid obstacles and the fluid region (hereafter 
known as the solid-fluid interface), n is the outward unit normal to u{x,t) is the 

velocity of the fluid flow, and 7 ^ 0 is the constant particle-adsorption coefficient. There is 
no adsorption when 7 = 0 , and the adsorption is instantaneous in the limit as 7 —)• 00 . 

Additionally, we assume that we have an incompressible Newtonian fluid, which satishes 
Stokes equations. 


—Vp + fj.V'^u=0, x£^f, (2.2a) 

^•11 = 0, x£^}f, (2.2b) 

u=0, X £ dQf, ( 2 . 2 c) 

where /r is the fluid viscosity and p{x,i) is the fluid pressure. 

2.1 Dimensionless equations 

We scale the variables via x = lx, u = Uu, i = {l/U)t, c = CooC, and p = {piLi/{5‘^l))p, 
where U \s a, characteristic velocity scale and Cqo is a characteristic concentration scale. The 
pressure scaling is chosen to balance the pressure gradient over the macroscale with viscous 
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Figure 1: Model schematic in two dimensions. Left: An example in which the macroscale 
porosity decreases in the direction of the flow. Right: A magnified view of a given cell uj{x), 
with microscale coordinate y G . 


forces over the obstacle microscale. Using these scalings in (12.ip . we obtain the dimensionless 
solute-transport equation 

dc 

— = V • (Pe“^Vc — uc) , X E ^f, (2.3a) 

—(5A:c = n • (Pe“^Vc — uc) , xEdi}f, (2.3b) 

where Pe = lAl/D^ the Peclet number, and k = 'y/{5Ll) are both assumed to be 0(1) 
parameters. As we see in ^ this choice yields a distinguished limit in which all mechanisms 
are present at leading order in the macroscale. We note that k may be expressed in terms of 
the Peclet number and the Damkohler number Da (which relates adsorption to diffusion) via 
k = Da/((5Pe). The parameter 5 that premultiplies the adsorption term in (|2.3bp is required 
to ensure a dominant balance between particle transport and adsorption, since an 0{6) 
adsorption rate over an 0(1/6) number of obstacles will lead to an overall 0(1) adsorption 
rate. It is shown in that an asymptotic analysis of a similar advection-reaction-diffusion 
equation with Pe = 0(1) and k = 0(1/6) leads to a breakdown in the asymptotic series 
considered. In the limit of large k, the solute is adsorbed very quickly and, for this work, 
could be tracked via an initial boundary layer in time of 0(l/fc) and, if necessary, diffusive 
boundary layers in space of 0(1/Vk) (though we do not consider these sub-limits further). 

The dimensionless version of ()2.2I) is 

—Vp + 6‘^V‘^u=0, xE^f, (2.4a) 

V • u = 0, xEOlf, (2.4b) 

u=0, X E dOf. (2-4c) 

In dimensionless units, the obstacles now form a d-cubic lattice of d-balls a distance of 6 
apart, and a d-ball with centre at x has radius 6R(x). A schematic of the dimensionless 
geometry is shown in the left-hand side of figured) 
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3 Homogenization 


In order to reduce the complexity of the problem geometry we homogenize the governing 
equations (IO])-(I22]) using the method of multiple scales. This will allow us to obtain 
effective equations on a simpler macroscale domain while formally capturing the relevant 
information about the microscale geometry. As is standard in homogenization theory, we 
introduce a microscale variable y = xjb and treat x and y as independent. The extra degree 
of freedom this allows is removed by imposing that the solution is exactly periodic in y, and 
hence any small variation between unit cells is captured through the macroscale variable x. 
As shown in figure [H the microscale variable y is defined in the unit cell uj{x), whereas the 
macroscale variable x spans the entire filter. The solid portion of the cell, occupied by the 
obstacle, is denoted by Us{x) and the fluid portion is 0 Jf{x) = uj{x)\uJs{x). The boundary 
of the unit cell determined by the obstacle is denoted by dujs{x), while the outer boundary 
of the unit cell is dijj{x). Further, we consider each dependent variable as a function of both 
X and y. Using the new variable y, we transform spatial derivatives in the following manner 

V I—)• Va, + yi (3-1) 

where Va, and Vy refer to the nabla operator in the x- and y-coordinate systems respectively. 

Our goal is to systematically reduce the geometrical complexity of the problem by deriving 
effective governing equations valid in the macroscale domain. For this purpose, the variables 
we use in our macroscale equations are quantities averaged over an entire cell lo(x}. To this 
end, we define the porosity 4>(x) to be 

^ K(a:)l. (3.2) 

since the unit cell has unit volume (|a;(a:)| = 1), and the volumetric average concentration C 
and volumetric average fluid velocity U (also known as the Darcy velocity in porous-media 
formulations) as follows 

C{x,t) = ] .. ( c{x,y,t)dy= I c{x,y,t)dy, (3.3a) 

l^(®)l J Lu(x) Jujf(x) 

U{x,t) = j u{x,y,t)dy= u{x,y,t)dy, (3.3b) 

l^(®)l j Ul(x) Ju)f{x) 

imposing that c = w = 0 in ujs{x). We note that, alternatively, we could have averaged 
over the fluid portion of the cell (using \ujf(x)\ in the denominator of (I3.3|) i rather than over 
the whole cell; this average, known as the intrinsic average, is most commonly used in the 
method of volume averaging . This would give the particle concentration in terms of the 
available fluid space of the membrane, rather than the overall particle concentration within a 
membrane. As we will see later, the volume average concentration is most convenient in our 
case since it allows us to directly deduce the effect of macroscopic changes in the porosity. 

We first consider the flow problem (|2.4p , since it has no dependence on the concentration 
field, and use the ensuing results in the solute-transport problem (|2.3p . 

3.1 Flow problem 

Under the spatial transform (|3.1I) . the flow equations (|2.4p become 


(d ^Vy -|- Va;) P -|- (Vy -|- dVa;) U — 0 , 

y G ujf{x), 

(3.4a) 

(Vy -|- 5Va;) • W = 0, 

y G ujf{x), 

(3.4b) 

u = 0, 

y G du}s{x), 

(3.4c) 

u, p periodic, 

y G du}{x), 

(3.4d) 
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where dojg is the obstacle boundary, and doj is the outer cell boundary as shown in figured) 
Expanding the flow velocity and pressure in powers of 5 as follows 


G{x,y,t) = Go{x,y,t) + 6Gi{x,y,t) + ... as 6 ^ 0, 

(3.5) 

where G G {u,p}, the flow equations ()3.4p at 0{6 ^) 

are 


-VyPo = 0. 


(3.6) 

Therefore, the leading-order pressure is independent of the microscale, i.e. po 
direct consequence of our pressure scaling choice in ^2.11 

The 0(1) terms in (13.41) are given by 

O 

II 

-VyPl -b Vylio = Va:P0, 

y G ujf{x), 

(3.7a) 

o' 

II 

O 

> 

y G ujf{x), 

(3.7b) 

Uo = 0, 

y G dujs{x), 

(3.7c) 

uq, Pi periodic. 

y G duj{x). 

(3.7d) 

Since (|3.7p is a linear problem, its solution can be written as 


Uq = -K{x,y,t)\/xPo, 


(3.8a) 

Pi = -n(a;, y, t) • V^Po + p{x, t), 

(3.8b) 

where the matrix function K and the vector function 11 satisfy the problem 


Id — Vyll -|- VyK — 0, 

y e Xf{x), 

(3.9a) 

II 

> 

y G Xf{x), 

(3.9b) 

K = 0, 

y G dujs{x), 

(3.9c) 

K, n periodic, 

y G dLjj{x). 

(3.9d) 


Here, Id is the d-dimensional identity matrix, and the dependence of K and 11 on a; is due 
to the dependence of the microscale boundary dujs{x) on the macroscale variable. 

To obtain the effective flow equation for the Darcy velocity (|3.3bp . we integrate (|3.8b i 
over ujf{x) to deduce 

U{x,t)= f u{x,y,t)dy ^ -K{^)Va,p, (3.10a) 

Jujf{x) 

at leading order, where is a scalar function defined by 

K{cj))ld= [ Kdy. (3.10b) 

We note that the integral of K in ()3.10bh is a multiple of the identity matrix due to the 
symmetry of the cell problem described by (13.Op . This would no longer be the case if, instead 
of a spherical obstacle, we had considered an anisotropic obstacle. Integrating (|3.4bll over ojf, 
and using the transport theorem (outlined in Appendix 1X1) in conjunction with the periodic 
and no-slip boundary conditions ()3.4cD - ()3.4dp . we obtain the following incompressibility 
condition for the macroscale flow velocity 

= 0. (3.10c) 

Therefore, the homogenization procedure of the Stokes velocity in the cell problem (12.41) via 
multiple scales yields Darcy’s law (I3.inp for the macroscopic velocity U as expected [27|. 
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3.2 Solute-transport problem 


We now perform a similar homogenization for the solute-transport problem. While the 
Dirichlet boundary conditions for the flow problem (j2.4b l are invariant to the multiple-scales 
transformation, this is not true for the transformation of the Robin boundary condition ()2.3b l 
in the solute-transport problem. Following 1^, 22], we address this issue by introducing the 
relation 


X{x,y) = R{x) -\\y\\, (3.11) 

where x(®,y) = 0 defines the solid-fluid interface^ Using the gradient transform (13.11) . the 
normal vector n in (|2.3bl) becomes 


Uy -(- xR 
riy -\- (5Va;i?|| 


(3.12) 


where riy = —y/||y|| is the outward unit normal on the obstacle boundary dujs{x), and 
6VxR accounts for the macroscale effect of varying obstacle size (see [l^ for details). 

Under the spatial transforms (13.11) and (13.121) . the solute-transport problem (12.3p becomes 


= i^y + <^V,) • (Pe-^ (Vy + <5V,) c - Sue) , 

-6'^kc = {uy + 6VxR) • (Pe"^ (Vy + JV^,) c - 6uc) + 0((5^), 
c periodic, 


y G 

(3.13a) 

y G dujs{x), 

(3.13b) 

y G du{x). 

(3.13c) 


Expanding c{x, y, t) in powers of 6 as in (|3.5I) and equating powers of 6, the 0(1) terms are 


0 = VyCo, 

y G u}f{x), 

(3.14a) 

0 - Tty * VyCQ, 

y G dujsix), 

(3.14b) 

Co periodic. 

y G duj{x). 

(3.14c) 

Therefore, we deduce that cq is independent of the microscale, i.e. cq = co(x,t). 
terms in (13.131) are 

The 0((5) 

0 = VyCi, 

y G w/(a:), 

(3.15a) 

^ x^O — 

, y G dujs{x), 

(3.15b) 

Cl periodic. 

y G duj{x), 

(3.15c) 

usinff the continuitv eauation for the flow (I3.7bl) 
allows us to write ci as 

. The svstem (13.151) is linear in 

Cl, which 

ci(a;, y, t) = -r(a;, y) ■ VxCo{x, t), 

(3.16) 

where the components of the function T satisfy 



o 

II 

<1 

y G xf{x), 

(3.17a) 

Tly^i - Tty • 

y G dusix), 

(3.17b) 

Fj periodic. 

y G du}{x), 

(3.17c) 


^Writing R{x) in (13.1111 assumes that it is a continuous function of the macroscale x, rather than a 
piecewise constant function evaluated at the centre of the relevant unit cell. This simplihes the subsequent 
analysis while affecting only the boundary condition (I3.13bll at higher orders than we need to consider. As a 
result, our final leading-order macroscale equation (I3.22all is unchanged by employing this simpliheation. 
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and Uy^i is the ith component of the unit vector Uy. 

The 0{5‘^) terms in (j3.13p are 

dc 

— = Vy • (Pe“^ (VyC2 + Va;Ci) - UiCq - UqCi) 

T Va; • (Pe (VyCi T VajCg) UqCq^ , 

-kco = Uy • (Pe“^ (VyC 2 + Va;Ci) - rriCo - uqCi) 

+ Va;-R • (Pe"^ (VyCi + VojCo) - uqCo) , 

C 2 periodic, 

Integrating (|3.18a|l over uj and applying the boundary conditions (|3.18bl) - (|3.18cp gives 


G ujf{x), 

(3.18a) 

G dujs{x), 

(3.18b) 

G doj{x). 

(3.18c) 

(IrMl-dT 

l 8 cP eives 


f ^ dy = [ Va, • (Pe ^ (VyCi + VajCg) - Uoco) dy 

JLJf{x) Jujfix) 


lu}f{x) 


Va,i? 


' duis{x) 


(Pe ^ (VyCi + VajCo) - lioco) ds - / fecods 

J du)a{x) 


(3.19) 


where ds denotes the differential element of the obstacle surface dujs{x). The first two terms 
on the right-hand side of (I3.19P can be simplihed using the transport theorem as stated in 
(HU) , giving 


lujf[x) 


dco 


Of{x) 


(Pe ^ (VyCi Va;Co) - ■Ugco) dy - f kco ds. (3.20) 

J dLda(x) 


Using that cg = cg(£c,t), \ojf{x)\ = (i){x) and (|3.16p . (|3.2np reduces to 

■ 1 / (Id-jf)dy Va,Cg - / Ugdy Cg ) - |5ws((/>)|A:cg, (3.21) 

Ot \ J 

where (Jr)ij = dVj/dyi is the transpose of the Jacobian matrix of P, the solution to the 
cell problem (|3.17p . From now on, we simply write 4> but it should be understood that the 
porosity will be, in general, a function of x. In a similar manner, we use ujf{4>) instead of 
ujf{x) henceforth. 

To express (I3.2ip in terms of the volumetric average concentration C{x,t) defined in 
(j3.3ap . we note that C{x,t) ~ cgdy = \ujf{x)\co{x,t) at leading order in 6. Using this 

relation and ()3.10ap . we find that 

Q(j f _ n _ \ 

^ = V, • iD{c/))V^C - - (U{ct>) + D{ct>)V^4>) \ - f{cf)C, (3.22a) 

at leading order in J, where the effective diffusion coefficient is 


D(</.)Id = Pe-M Id - - / J?dy , 


and the effective adsorption coefficient is 




<(> 


1 - 


(3.22b) 


(3.22c) 


using the fact that <^(») = 1 — VdR{x)'^ for spherical obstacles, where is the volume of 
the unit ball with dimension d, so V 2 = and V 3 = 47 r/ 3 , and \du}s{4>)\ = dVdR{xY~^. 
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(f) ^ _ 

(a) Log-plot of the permeability K{(j)). (b) The diffusion coefficient D(4>). 


Figure 2: The functions (a) and (b) PeD{(j)) (defined in (jli.inbp and (j3.22bp respec¬ 

tively) calculated using Comsol Multiphysics, for ball obstacles whose centres are located in 
a cubic lattice, but whose radii can vary as a function of space. Due to the cubic lattice 
arrangement, the minimum porosity is 1 — 7r/4 ~ 0.21 in two dimensions, and 1 — tt/G ~ 0.48 
in three dimensions. Note that K ^ oo as 4> ^ 1, but we bound the graph for illustrative 
purposes. 


We recall that the Darcy velocity that appears in (j3.22ap is given by U{x,t) = —K{4>)Vxp, 
with the effective permeability K defined in (j3.10bl) . As for the permeability K, we note the 
right-hand side of ()3.22bh is a multiple of the identity matrix due to the symmetry of the 
cell problem (|3.17l) . 

The effective permeability K{(f>) and the effective diffusion D{(p) can be computed by 
solving the respective cell problems (j3.9p and (I3.17P for a given cell porosity cj), determined 
by the size of the spherical obstacle. We do this numerically using the finite-element software 
Comsol Multiphysics and find that both K and D increase as the porosity increases, as 
expected (see figure [2]). The additional degree of freedom afforded by the three-dimensional 
problem means that the permeability and effective diffusion are larger in three than in two 
dimensions. Our calculations show that the leading-order effective diffusion coefficient is not 
affected by the fluid flow, hence the effective diffusion shown in figure [2](b) is equivalent to 
that calculated in [^, which considered the limit oi U —)• 0, A: —)• 0. 

The homogenization procedure has significantly decreased the geometrical complexity of 
the multiply-connected domain in the full problem (|2.3p . while only slightly increasing the 
complexity of the governing equations for the homogenized problem (|3.22p . This complexity 
manifests in the variable coefficients of the governing equation; these coefficients reflect the 
microscale structure of the problem in a systematic manner. 

In addition to the increase in computational efficiency, our homogenized model gives 
insight into the physical effect of a non-uniform porosity. Namely, the homogenized governing 
equation (|3.22all indicates that a gradient in the porosity enters the equation as a term similar 
to the advection term due to the Darcy velocity. In other words, a non-uniform porosity in 
the membrane gives us an extra degree of freedom with which to either enhance or reduce 
the effect of a fluid flow on the solute concentration distribution. In the next section we 
determine how exactly to vary a macroscale porosity gradient to improve the effectiveness 
of a membrane filter. 
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4 A uni-directionally graded filter 


4.1 Model set-up 

In this section we use the homogenized equations to understand and quantify the effect of 
porosity gradients on filter efficiency. We model a filtration process using (|,'I.22I) in the steady 
state, i.e. setting dC/dt = 0, coupled with appropriate boundary conditions at the inlet and 
outlet that represent two reservoirs. Since we are interested in the interaction between the 
fluid flow and the porosity gradient in the advection term, we consider a filter that is graded 
in the same direction as the fluid flow. Whilst it can be difficult to fine-tune filter porosity 
due to manufacturing constraints, it is acknowledged that a negative porosity gradient can 
improve filter efficiency. As discussed in the Introduction, a quantitative understanding of 
this behaviour is an open question. 

Specifically, we consider an industrial process whereby a filter with a uni-directional 
constant porosity gradient separates two reservoirs, and a flow is induced through the filter 
in the same direction as the gradient in porosity. We are interested in the steady-state 
operation of this process, and how the porosity gradient affects particle adsorption. We 
define the direction of porosity gradient as x, where x € (0,1) within the filter (since the 
dimensional characteristic length I is chosen to be the length of the hlter). Thus the set-up 
is similar to that illustrated in figure [TJ The upstream is defined for x G (—oo, 0) and the 
downstream for x G (l,oo). We emphasize that the full problem occurs within a three- 
dimensional domain (and take d = 3 henceforth) and it is only the variation in porosity that 
is one-dimensional. 

Provided the boundary conditions allow for uni-directionality, the Darcy equations (I.S.lOp 
yield a uni-directional constant macroscale flux U = Ue^ (where is the unit vector in 
the x-direction), at all points in the filter, regardless of the porosity. We are able to take 
U = 1 without loss of generality by choosing the characteristic velocity scale U to be the 
average velocity in the x-direction, a quantity that is easily measured experimentally. The 
homogenized governing equation for the concentration C within the filter (I3.22ap is then 


dx 


,,dC C 




dx (p 




dx 


Considering the limit of (|3.22bp and (14.11) as p — 
upstream C~ and downstream concentrations, 


dx 

dx 


Pe 


Pe 


-1 


dC- 
dx 
_^dC+ 
dx 


- c-j = 0, 

- cA = 0, 


fiP)C, XG(0,1). (4.1) 

provides governing equations for the 


X G (—oo, 0), 

(4.2a) 

X G (1, oo). 

(4.2b) 


We impose continuity of concentration and concentration flux at the boundaries between 
the filter and the reservoirs. In the far-field of the reservoirs, the concentration tends to 
a constant value. We may take the upstream concentration C~ —>■ 1 (by choice of our 
nondimensionalization), while the downstream concentration C~^ tends to a constant value 
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that must be determined as part of the solution. Mathematically, this corresponds to 


c- 1, 


Pe 


idC+ 


dx 


-c..,,.,a). 


dC+ 

dx 


0 , 


Equations (14.21) . (I4.3al) and (|4.3fl) give 

C~ = 1 — Aexp(Pex), C~^ = B, 


X —oo, 

X = 0, 

X = 0, 

X = 1, 

X = 1, 

X oo. 


(4.3a) 

(4.3b) 

(4.3c) 

(4.3d) 

(4.3e) 

(4.3f) 

(4.4) 


where A,Bg [0,1] are two constants to be determined. These expressions may be substituted 
into the remaining boundary conditions (l4.3bP " (l4.3ep to yield two boundary conditions for 
C, thus closing the elliptic problem (14.11) . 


D{cP) 


dx 


C 

4 > 


(l + D(0)^) --1. 

^ _ Q 

dx I?!) dx ’ 


X = 0, 

(4.5a) 

X = 1. 

(4.5b) 


We solve the system for the filter (14.11) and (|4.5I) numerically using a second-order accurate 
finite-difference scheme with 10^ grid points. To ensure that our results are experimen¬ 
tally relevant, we make comparisons between simulations that use the same transmembrane 
pressure difference. This corresponds to modifying the dimensionless parameters Pe and k 
between simulations, and is discussed in more detail in Appendix [Bj 
We consider a linear porosity model of the form 


4>{x) = (j)o + m{x — 0.5), (4.6) 

for various mean porosities (po and gradients m. This corresponds to a uniform filter when 
m = 0. As we expect, the concentration C decreases with x as the contaminants are trapped 
by the filter medium (figure [3^). We also observe that the rate at which the concentration 
changes in x (which is related to the rate of particle removal) decreases as the porosity 
gradient increases (figure [3^). We probe this result further in the next section, to explore 
which filter set-up is most efficient at contaminant removal. 


4.2 Filter efficiency 

Since the goal of filtration is to maximize particle adsorption, it is helpful to consider the 
particle adsorption distribution, described by the right-hand side of (BH), within the filter. 
To this end, we introduce 


A(x) =/(())(x))C(x). (4.7) 

This measures the concentration of contaminant removed per unit time at position x in 
the filter. For the concentration profiles depicted in figure [3^, we show the corresponding 
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(a) Contaminant concentration distribution. (b) Adsorption distribution. 


Figure 3: Numerically determined concentration and particle removal profiles for the porosity 
function (j4.6p using = 0.75 with three different gradients, m = —0.3,0,0.3. We use the 
reference valnes (/>o = 0.75, m = 0, Pe = 3, /c = 1 from which to modify appropriate 
parameters (as discnssed in Appendix [B]) . The concentration C{x) is obtained from (14.ip 
and (14.5p using a second-order accnrate finite-difference scheme, from which the adsorption 
distribution A(x) can be deduced using ()4.7I) . 


A in figure [HJa for comparison. A natnral measnre to evalnate filter efficiency is the total 
contaminant removal by the filter per nnit filter area per nnit time, given by 


T = 



Using (14. ip . ()4.3dp . and (Ii31) we see that 


(4.8) 


r = l-C'+<l, (4.9) 

and thus total contaminant removal is directly related to outlet concentration, as expected. 

We investigate how a porosity variation affects the total particle removal within a filter by 
calculating T as we vary the mean porosity (/>o and gradient m in the linear porosity function 
(IMI). Although T is observed to increase with decreasing (j)Q, as we would expect, there 
are diminishing returns as we continue to decrease (j)Q (see figure 0^). However, we observe 
only a weak dependence on m (figure Ht'-) whereby, upon including a porosity gradient, there 
is an increase in the total particle adsorption T. Interestingly, we find that the sign of m 
does not affect T, exhibited through the symmetry about m = 0. Whilst it may appear 
surprising that the orientation of the filter has no effect on total particle adsorption, this 
phenomenon arises as a result of the commutativity of the obstacle layout: the concentration 
after each subsequent obstacle is a proportion of the incoming concentration. Hence, the 
contaminant removal by the entire filter is a product of all these proportions and, moreover, 
any permutation of the sinks will result in the same total contaminant removal. If our model 
included pore blocking as a result of particle adsorption, this symmetry feature would no 
longer hold. 

Thus, while the metric T is useful, it cannot account for the superior performance often 
observed for filters with a negative porosity gradient. To explain this experimental observa¬ 
tion, we must investigate A further. Although long-term blockage is inevitable, a particular 
aim of an efficient filter is to avoid local blockages while the rest of the filter is still functional 
i.e. to have the filter fail ‘all at once’ rather than in one place. Thus, while the total con¬ 
taminant removal is important, we must also consider the uniformity in which it is removed 
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(a) T with modified parameter values. 



(b) S with modified parameter values. 


Figure 4: Plots of the filter efficiency metrics T and S, defined in (j4.8l) and (|4.1U|) respectively, 
for varying (po and m in (p{x) = (po + m{x — 0.5). Each grey curve represents a different value 
of 4>q, which increments in steps of 0.05 from (f>Q = 0.65 to (j)Q = 0.9. We vary m for a given 
4>o such that 4> € [0.55,0.95] (which falls within the allowable range of </> G [1 — vr/6,1] in 
three dimensions). Therefore, the available range of m varies with (pQ. We use the reference 
values (j) = 0.75, Pe = 3, /c = 1 from which to modify appropriate parameters (as discussed 
in Appendix [B]) . The asterisks in (b) denote the lowest value of S for a varying m with a 
fixed value of (pQ. The dashed black curves denote the asymptotic approximations of T and 
S from results derived in ^4.31 and offer outstanding agreement with the numerical results 
at a fraction of the computational cost. 


by the filter. To quantify this, we introduce the metric 

S = [ |A(x) -Tj dx; (4.10) 

Jo 

S is minimized (and zero) for uniform contaminant removal and so we posit that a smaller 
value of S corresponds to a more efficient filter. 

In contrast to our analysis of the total particle adsorption, T, we find that S strongly 
depends on the orientation of the filter (see figure [Da). Additionally, a negative porosity 
gradient m yields a concentration profile that is closer to providing spatially uniform con¬ 
taminant removal than a positive m. We note however that the minimum value of S for a 
given mean porosity (po (shown as asterisks in figure (D?) does not necessarily correspond to 
choosing the porosity gradient as negative as possible. This feature occurs because a porosity 
gradient that is too negative will have less particle adsorption near the entrance to the filter 
than towards the exit, thus also creating an uneven adsorption distribution. 

Our analysis therefore indicates that, for a given mean porosity, the porosity gradient 
has a much larger effect on the distribution of particle adsorption within the filter than it 
does on the total particle adsorption. This provides a quantitative interpretation for the 
experimental observation that efficiency in contaminant removal is improved when a filter 
with a negative porosity gradient is used. As the shape of the concentration profile relates 
to the problem of localized blocking, we hypothesize that a negative porosity gradient allows 
for a more uniform contaminant removal. This increases the time until localized blocking 
becomes a problem, and thus increases filter efficiency. 
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4.3 Slowly varying porosity 

The results we have presented so far provide a very useful quantitative analysis on the role 
of porosity gradients in filtration. However, because they were obtained numerically, it is 
still rather cumbersome to perform major parameter sweeps through the space of porosity 
functions and the dimensionless parameters Pe and k (although we emphasize that results 
we presented in the previous section are still orders of magnitude cheaper than solving the 
full equations in a complex geometry). 

In this section, we take a further simplifying step by exploiting the fact that manufactur¬ 
ing constraints limit the possible variations in porosity. We determine approximate analytic 
solutions to the equations presented, significantly increasing the computational efficiency and 
allowing us to gain insight into the concentration profiles we obtain. To implement this, we 
explore the limit of a weakly varying porosity such that d(p/dx = 0(e), where e <C 1. 

We expand C and ()>(x) in powers of e as follows 

C'(x) ~ C'o(x) + eC'i(x), (/)(x) ~ ())o + as e ^ 0, (4-11) 


where (po is the mean porosity, thus constant in x. We expand D(4i{x)) and f{4>(x)) in the 
same manner as (p(x), i.e. where the leading-order term is constant. The coefficients are 
given by 


Do — D((f)o), 

d(l-<f>o)( Vd 


h = k- 


4>o 


( Vd y/" 
Vi-</'oy ’ 


Di(x) = (j)i{x)^((t)o), 
ji(x) = k(t)i(x)—j- 


Vi-</>oy 


(4.12) 


As the porosity function is an input to our model, the asymptotic expansions for (l)(x), 
D((j)(x)), and /(())(x)) are known. For the linear porosity functions (|4.6I) considered in the 
previous section, we have e = m, ())i(x) = x — 0.5, and the asymptotic limit we consider 
is equivalent to m —)• 0, i.e. close to uniform porosity. Note that the theoretical maximum 
value of e for a linear porosity function is vr/O ~ 0.52 in three dimensions. However, the 
maximum value of e will be smaller in practise. 

Substituting the asymptotic series ()4.11l) into the system (j4.ip . ()4.5p and equating powers 
of e, the 0(1) system is given by 


d (^dCo Co 


dOo Co 


dx \ dx 4>o 
_ C 

dx 


dx (j)o 


The system (j4.13l] is solved by 


= foCo, 

X G (0,1), 

(4.13a) 

= -l, 

X = 0, 

(4.13b) 

= 0, 

X = 1. 

(4.13c) 

ab{x — 1) 

— sinha5(x — 1)] , 

(4.14a) 


where 


1 


1 


+ (i + 42),i„h„6 + 26cosh<.(,' 


(4.14b) 
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Figure 5: The volumetric concentration profiles generated from numerical (solving the sys¬ 
tem (I4.ip . (I4.5p using a finite-difference scheme) and asymptotic (obtained from (|4.14p 
and (j4.16p i methods. The parameter values used are Pe = 5, A: = 1, 4>{x) = 0.75-|-m(x —0.5). 
The two-term asymptotic solution is indistinguishable from the numerical solution in (a). 


dx 


The 0(e) terms in the system (j4.ip and (14.5p are given by 

dOi Cl 


On 


O, 


dx (po 

dCi Cl 

^ dx (po 


(4.15a) 

dOi Co dcPi 


dx (po dx ’ 

Using variation of parameters, the system (|4.15p is solved by 


x = 0, (4.15b) 
X = 1. (4.15c) 


^ exp(ax) 
^ sinha 6 


(Ai(x) -|- aBi) sinhaftx -|- {A 2 {x) + 01 B 2 ) sinha 6 (x — 1) 


Here, 


1 /■ 

yli(x) = —/ ^(t, Co(t)) exp(—at) sinh a 6 (t — 1 ) dt, 
abDo Jx 


1 O 

A2{x) = —/ ^(t, C'o(t)) exp(—ot) sinhaftt dt, 
abDo Jo 

W b(^A2{l)-2^)_' 

J VMi(0) + iV(o,Oo(0))/(aOo), 


Oi 

O 2 


abDo 

—q b 
b -q. 


d 


g{x, C{x)) = C{x)) -h /i(x)C'(x), 

dO C'(x) (— d(pi (pi{x) 


N(x,Ca)) = D,(x)—- 

q = sinh 06 -|- 6 cosh ab. 


Do 


dx 


4>o 


(4.16) 

(4.17a) 

(4.17b) 

(4.17c) 

(4.17d) 

(4.17e) 


(4.17f) 

In addition to the increase in computational efficiency, this analytic result also allows us 
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to determine closed-form approximations up to 0(e) for the metrics T and S, 

5 ~ / |/o^^o(x) + e (Mx)Co(x) + /oCi(x)) - T| dx, (4.18b) 

Jo 

The asymptotic concentration profile (up to two terms) agrees very well with the profile 
obtained numerically, even when the porosity gradient takes physically extreme values (see 
figure E]). Further, the closed-form nature of (|4.18p allows us to produce approximate results 
for T and S (shown as dashed black lines in figure U]) at a fraction of the computational 
cost. As with all the analytic results in this section, we observe excellent agreement with the 
numerical results. 


5 Discussion 


We have systematically derived a macroscopic model for a porosity-graded filter from mi¬ 
croscale information by generalizing standard homogenization theory for near-periodic sys¬ 
tems. The result is an advection-diffusion-reaction equation for the solute concentration 
within the filter as a function of the porosity distribution and operating conditions. In 
particular, this equation has allowed us to investigate how porosity gradients affect solute 
trapping and filter efficiency. To this end, we have defined two suitable metrics to quantify 
the implicit effect of blocking. The first corresponds to the rate of total contaminant removal. 
The second measures how close to uniform the adsorption of the contaminant is within the 
filter and is an indication of the propensity of a filter to localized blocking. By performing 
a computationally efficient parameter sweep, we have been able to determine and quantify 
that, for a given mean porosity, a porosity gradient has a much larger effect on the second 
metric than the first metric. In general, we have found that a filter with a decreasing porosity 
will be less prone to localized blocking and we have quantified the optimal gradient as a func¬ 
tion of the operating conditions. This allows us to understand the experimental observation 
that a decreasing porosity can lead to a greater filter efficiency; the porosity gradient lowers 
the risk of localized blocking, while maintaining the rate of total contaminant removal. 

The homogenization procedure we used allowed a near-periodic, rather than a strictly 
periodic, microstructure to be considered. The macroscale variation that near-periodicity 
allows is a vital feature of a porosity-graded filter. The resulting macroscale equations are 
computationally cheap to solve, allowing us to efficiently explore the experimental parameter 
space to quantify filter efficiency. 

With regards to an ‘ideal’ filter, we note that it is impossible to globally optimize both 
of the two metrics (T and S) that we created. Indeed, we found that porosity functions that 
increase total adsorption were further from uniform adsorption. Thus, there is a trade-off to 
be made between maximizing adsorption and minimizing potential blockage issues, an idea 
supported in the network model for trapping 2^. Whilst the question of determining an 


‘optimum’ porosity function will depend on the specifications of the end-user, our work pro¬ 
vides an inexpensive method to increase filter efficiency given these requirements. However, 
our model does suggest that, in the absence of additional constraints, it is wiser to use a 
filter with a negative rather than positive porosity gradient. While this will not affect the 
total adsorption, it will provide a more uniform adsorption rate throughout the filter, which 
should reduce blocking issues. 

This work can also be used to investigate filters comprised of a series of material layers, 
each with a different mean porosity. One way to achieve this would be to derive suitable 
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boundary conditions to couple each of the regions, using the method in ^4.11 Another 
(simpler) idea, would be to use a differentiable approximation of the porosity {e.g. cubic 
splines defined between the centres of neighbouring layers) in the macroscale model defined 
in the first part of this paper. 

Although we considered a filter whose microscale structure consisted of balls, it is possible 
to consider an arbitrary shape for the microscale obstacle. However, we chose our partic¬ 
ular microstructure to allow an explicit macroscale equation, and thus significant analytic 
progress to be made, and the same will not hold for the general microstructure case. It is 
a simple task to combine this work with that of [l^, who considered a general curvilinear 
coordinate transform to map a near-periodic microscale to a periodic domain, thus allowing a 
homogenization method to be applied. However, as the coefficients within the resulting gov¬ 
erning equations would have to be determined at each point in space, the ‘reduced’ equations 
are still computationally expensive. 

One drawback of the method presented here is the restriction to a near-periodic mi¬ 
croscale, meaning that materials that vary wildly on the microscale cannot be treated sim¬ 
ilarly. A viable extension to this work would be to test how robust the assumptions on a 
near-periodic microscale structure are, and how far away one could get before the macro¬ 
scopic description we derived broke down. Although the near-periodic assumption is required 
to derive the macroscale equations, it may be the case that it is only the average properties 
of the microstructure that must possess this property. This could be tested by performing 
numerical simulations on the full problem for randomly placed spheres. Similar simulations 
have been carried out in 2^ for diffusing particles only, and it was found that the concen¬ 
tration distribution for a random and cubic array of spheres agree in the high-porosity limit. 
On physical grounds, we expect the same result for our work in the high-porosity limit. 
Furthermore, in the same limit, we expect the results for a general microscale obstacle to 
coincide with the high-porosity limit of the governing equations derived in this paper, using 
appropriate results for pore volumes and obstacle surface areas. 

One feature that is lacking in our model is the dynamic effect of blockage (membrane 
fouling). As blocking would temporally change the geometry of the problem, the flow and 
particle transport problems would be fully coupled, and hence difficult to investigate. Despite 
the lack of explicit blocking in our model, we were able to implicitly infer the long-term 
effect this would have by quantifying the particle adsorption distribution within the filter. If 
blocking was also considered, we would expect the results in this paper to provide a quasi¬ 
static description of that model. Additionally, the linear trapping rate that we imposed is a 
generalized approximation of a feature that will depend on the solute solubility and the hlter 
material adsorption. Intuitively, one would expect a Michaelis-Menten type dependence for 
the adsorption rate on the number of particles, i.e. approximately linear for a small number 
of particles and bounded above for a large number of particles. Our approach is valid because 
of the dilute suspension assumption but, if one is in the same asymptotic limit, it is a trivial 
task to extend the analysis in this paper to a general nonlinear adsorption rate if required. 

Finally, we note that this work has the potential not only to guide filter manufacture and 
operating conditions, but also to provide assistance to other industries that use functionally 
graded materials, such as heterogeneous artificial body tissue in tissue engineering [^, [29| . 
or graded electrodes in lithium-ion batteries 30|. Furthermore, as the technology to produce 
functionally graded materials grows 3l|, the potential experimental parameter space will also 
increase. Forming appropriate mathematical models and maximizing the analytic progress 
that can be made will significantly expedite the exploration of this parameter space and 
result in faster technological growth. 
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A Transport theorem for a varying domain 


During the homogenization procedure, we use the transport theorem to evaluate integrals 
over microscale cells which vary in the macroscale variable x. To do this, we must determine 
the rate of change of the fluid domain ujf{x) with respect to x. 

We note that ojf{x) has a fixed outer boundary, doj, and an interior boundary on the 
surface of the obstacle of radius R{x), dus{x). The rate of change of dujs{x) with respect to 
X is VxR. To see this, consider the difference between integrals over the domains ujf{x + ^ei) 
and ojf{x) as ^ ^ 0. The resulting domain of integration is a shell whose thickness is ap¬ 
proximately ^dR/dxi as ^ 0. As we are considering an integral over u}f{x) = uj{x)\ujsix), 

the relevant velocity of the interior boundary is — 

Therefore, the transport theorem for the geometry considered in this paper is as follows. 
Let v{x,y,t) be a continuous vector field and periodic in y on the outer boundary duj{x). 
Then 

Vx- v{x,y,t)dy= Vx • v{x,y,t)dy - VxR • v{x,y,t) ds. (Al) 

Jijf{x) JduJs{x) 


B Parameter modification 


Whilst the two dimensionless parameters Pe and k are useful for a mathematical analysis of 
the governing equations, within an experimental set-up it is not these parameters that will, 
in general, be kept constant as we vary the membrane porosity. When analyzing results from 
our model, it is useful to consider mathematical variations that correspond to experimental 
variations in order to maximize the experimental relevance of our work. 

In practice, the control parameter within a filtration system is typically the transmem¬ 
brane pressure difference (note that the same calculation applies if the volumetric flux is the 
control parameter that is used instead). This corresponds to a linear variation in the average 
velocity in the x-direction, U. For a given (dimensional) transmembrane pressure differ¬ 
ence AP, we can solve (I3.10h l and (I3.10b l in one dimension to deduce that (in dimensional 
variables) 


U[(j){x)] 


{SlfAP / dx \ ^ 
y Vo K{4>{x))) 


(Bl) 


Here, the square brackets denote that is a functional of </>(x). 

Therefore, we are able to determine how U[(p] varies (relative to a reference value for 
the transmembrane pressure difference) for different porosity functions 4>{x). That is, for a 
given reference function (pj.ef{x), and parameter values Pe[(/)ref] and k[4>,-ei], we can determine 
the effect of changing 4>{x) on whilst keeping the transmembrane pressure difference 
constant. This, in turn, allows us to calculate the values of Pe[(/)] and k[(l)] in the dimensionless 
model that correspond to the new function (j) while keeping the transmembrane pressure 
difference constant. Accordingly, when we discuss S and T in ^4.21 we take Pe[(/)] and k[(j)] 
to vary appropriately, relative to a reference point, defined by 4>o = 0.75, m = 0, Pe = 3, 
k = 1 (illustrated in figure [6]). 


References 

[1] A. W. Zularisam, A. F. Ismail, and R. Salim. Behaviours of natural organic matter in membrane 
filtration for surface water treatment - a review. Desalination, 194(1):211-231, 2006. 


19 








10 


2.5 


Pe 


-0.4 


Increasing (po 


- 0.2 


0 

m 


0.2 


0.4 


(a) Variation of Pe. 


k 


2 



Increasing 4>o 

-8.4 -0.2 0 0.2 0.4 

m 


(b) Variation of k. 


Figure 6: Plots of how the dimensionless parameters Pe and k vary in a manner consistent 
with keeping the transmembrane pressure difference constant between experiments for vary¬ 
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